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Abstract 

The leading order vacuum polarization contribution of very light millicharged 
fermions and scalar (spin-0) particles with charge ee and mass /x to the Lamb 
shift of the hydrogen atom is shown to imply universal, i.e. ^-independent, upper 
bounds on e : e ;$ 10~ 4 for u ;$ 1 keV in the case of fermions, and for scalars 
this bound is increased by a factor of 2. This is in contrast to expectations based 
on the commonly used approximation to the Uehling potential relevant only for 
conventionally large fermion (and scalar) masses. 



The recent observation [1] of an optical rotation of linearly polarized laser light gener- 
ated in vacuum by a magnetic field may be due [2] to the photon initiated pair production 
of very light charged fermions with mass /i ~ 0.1 eV and charge ee where e ~ 10~ 6 . There 
exist, however, very strong astrophysical, cosmological and laboratory constraints [3, 4, 5] 
which exclude the quoted values of e and /i. Some of these constraints may nevertheless be 
relaxed in specific paraphoton scenarios [5]. Further constraints on the mass and charge 
of light charged fermions may be obtained from their leading order vacuum polarization 
contribution [3, 7, 8] to the Lamb shift of the hydrogen atom 



or from their higher order contribution [7] to the anomalous magnetic moment of the 
muon. The Lamb shift constraints for conventionally large fermion masses studied in 
[3, 7, 8] were, as ususal, based on the commonly used approximate leading order vacuum 
polarization contribution 



However, the upper bound presented in [3, 9] is not appropriate for fi <, 1 keV. Even 
much lower values of fi, fi < 1 eV, were considered in [2], for example, and it is therefore 
necessary to study the consequences of going beyond this standard approximation relevant 
only for fi > am e . 

The lowest order Coulomb interaction V = —a/r for a point nucleus is modified at 
the 1-loop level according to V = —a/r + 5V where [10, 11] 



and a = e 2 /4n = 1/137.036. This leading order Uehling potential yields, instead of the 
approximate equation (2), the exact leading order expression 



SE = E(2S 1/2 ) - E(2P 1/2 ) 



(1) 
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with the normalized radial hydrogen wave functions R n i given by 
where p = r/a and a -1 = ccm e . The integration over r in (4) yields 
where 

I (a*) = l°° du ( 1 + /^~!, (7) 

V ; A V 2 «V («* + 2«) 4 v ; 

with the effective coupling a* = am e / //. Requiring this leading order vacuum polarization 
contribution \5Ey P \/27rh to the Lamb shift not to exceed O.Of MHz, corresponding to a 
la discrepancy between the measured and calculated shifts [12, 13], yields the constraints 
on e shown in Fig. 1. (Imposing a la error, i.e. increasing the experimental uncertainty 
to 0.02 MHz, increases the upper bound on e by a factor of y/2 .) It is interesting to 
note that, contrary to what one might naively expect, it turns out that the integral I (a*) 
in (6) entirely compensates the large enhancement factor a* 2 already at p: <, 1 keV due 
a* 2 1 (a*) = for a* 3> 1. (This asymptotic result is easily obtained by observing that 
the maximum of the integrand in (7) moves to oo as a* — > oo where the integrand reduces 
to u(a* + 2w)~ 4 .) Therefore the upper bound in Fig. 1 saturates below /i ~ 1 keV which 
implies a rather moderate but universal, i.e. /^-independent, upper bound 

e < 1.085 x 10" 4 for /i < 1 keV . (8) 

Although in the different context of electronic vacuum polarization corrections to the en- 
ergy levels of exotic atoms, the exact result (7) and its asymptotic limit was previously 
presented [14] in a different form. 1 Despite the fact that our eqs. (6) and (7) are particu- 
larly suited for directly understanding the reasons for the universal bound on e discussed 



1 The integral in (7) as well as the ones for the (physically irrelevant) individual 2S and 2P contributions 
can in general be calculated analytically. The rather lengthy analytic expressions of the latter ones, 
together with their asymptotic limits (a* 3> 1), are given in Appendix A of [14] by identifying k,i — a* /2. 
Although these expressions are less transparent for understanding our universal (scaling) upper bound 
on e, it should be mentioned that the individual 25 and 2P contributions do not scale asymptotically 
but increase as In a* — ^ and In a* — |, respectively, in units of — a 3 e 2 m e /6n. The advantage of our more 
direct approach is that it can be immediately applied also to the case of scalar particles studied below. 
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above, we also present for completeness the analytic result for the integral in (7): 



with K = a* 12 and L = ln(K +f\ =T) relevant for k > 1, and L — for k < 1. This 



analytic result has originally been given explicitly in [15] in the context of mesonic atoms. 
For a comprehensive review on radiative QED corrections to hydrogenlike atoms we refer 
the reader to [16]. 

The dashed curve in Fig. 1 displays the upper bound as obtained from using [3, 8] 
the approximate expression (2) which, imposing \8Ey P \/2irh < 0.01 MHz, implies e < fx/ 
(26.62 MeV) as was assumed to hold for \x > 1 keV [8]. It can be already seen from Fig. 1 
that this approximate bound is not correct for \x much below 10 5 eV. To illustrate this 
more clearly, we confront in Fig. 2 the approximate and exact upper bounds on e, as 
obtained from (2) and (6), respectively, using a linear scale for e. Their convergence at /i 



It should be mentioned that corrections to the anomalous magnetic moments of elec- 
trons and muons, or to the hyperfme splitting, for example, are genuinely swMeading and 
thus the resulting bounds are much weaker [7]. Furthermore, considerations of the level 
shifts of exotic atoms [14] obviously can not improve the upper bounds on e as obtained 
for the hydrogen atom, since the associated relative theoretical uncertainties are generally 
more significant [17]. 

Next we also calculate the exact upper bound for light millicharged scalar (spin-0) 
particles of mass In this case the standard leading order fermionic QED con- 
tribution x(l — x) In [1 — x(l — x) q 2 /fi 2 ] appearing in the integrand of the Feynman- 
parameter integral in the vacuum polarization tensor [10, 18] has to be replaced by 





> 10 5 eV is due to I (a*) 



j_ 

40 



for a* -> in (7). 
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|(2x — l) 2 ln [1 — x(l — x) g 2 //i 2 ]. Performing now the exact integrations as outlined in 
[10], for example, one arrives at 



SV.(r) = - a -?fr *.«-<"« (l - A) (10) 

for the 1-loop correction of the leading order Coulomb interaction V = — ct/r + 5V^. The 
calculation of this vacuum polarization contribution to the Lamb shift is now analogous 
to the above case of fermions, and the exact leading order contribution of scalars reads 

5E^ P = -^e 2 a* 2 I s (a*) (11) 

where 

1 \ Vu^ : 



IJa*) = du 1 - — -r , (12) 

K J J l V J («* + 2w) 4 V ; 

to be compared with the fermionic result in (6) and (7). Imposing the same experimental 
upper bound on \5Ey P \/2nh of 0.01 MHz as above for fermions yields the constraints on 
e shown in Fig. 1 for scalars. The scalar upper bound in Fig. 1 again saturates below 
/j, ~ 1 keV due to a* 2 I s (a*) = ^ for a* 3> 1, which is the same asymptotic result as in 
the fermionic case for a* 2 1 (a*). Therefore, since the scalar contribution in (11) is 4 times 
smaller than the fermionic one in (6), the universal, i.e. /^-independent, upper bound for 
the charge of scalar particles is a factor of y/i larger than the one given in (8), i.e. 

£<2.17xl0~ 4 for /i<lkeV. (13) 

The dashed curve in Fig. 1 for scalars displays the upper bound as obtained from utilizing 
the approximate expression [7] 

5 / \ 2 
240V 

which is obtained from (11) using I s (a*) — > ^ for a* — > 0. (Notice that this result is 
a factor of 8 smaller than the fermionic approximate result in (2).) It can be seen from 
Fig. 1 that the approximate scalar bound is not correct for /j, much below 10 5 eV as in the 
fermionic case. Due to the different integrand in (12), as compared to (7), the difference 



SEb*-^ ^) e 2 (14) 



between the approximate and exact upper bounds is more pronounced than for fermions 
which is illustrated more clearly in Fig. 3 (to be compared with Fig. 2). Finally it should 
be mentioned that the scalar scenario is experimentally somewhat favored (cf. Table V of 
over the fermionic one. 



Again, as in the fermionic case, our eqs. (11) and (12) are particularly suited for directly 
understanding the reasons for our universal bound (13) on e. Nevertheless we also present 
for completeness the analytic result for I s in (12) as obtained from a straightforward 
integration: 



a * 2ls{a * ] = 8^ 



K ( 20 9 1 A 

2p L T)( 16 - 24K2 + fc4 ) i 



(15) 



with k and L as in (9). 



To conclude, we have shown that the contribution of very light fermions with charge ee 
and mass fi to the Lamb shift of the hydrogen atom implies a universal, i.e. /^-independent, 
upper bound e £ 10 -4 for fi <, 1 keV. This result is only obtainable by utilizing the 
exact expression for the Uehling potential rather than its standard approximation [11, 18] 
commonly used in QED. Since the scalar spin-0 scenario is experimentally somewhat 
favored [2], we calculated the exact upper bound for millicharged light scalar particles as 
well, which turns out to be a factor of 2 larger for /i <, 1 keV than for fermions. 
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Figure 1: Allowed upper bounds of e according to the exact leading order contribution 
(6) and (11) of fermions and scalars, respectively, to the Lamb shift of the hydrogen atom 
corresponding to a la discrepancy of 0.01 MHz between theory and experiment. The 
dashed 'fermion' curve corresponds to the bound suggested in [3, 8] being based on the 
approximate expression (2) as discussed in the text. The dashed 'scalar' curve is obtained 
from the approximate expression (14). (In order to avoid any confusion it should be 
stressed that the area above the respective curves is excluded.) 
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